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ABSTRACT 


Seasonal,  monthly,  weekly,  and  hourly  sea  level  variability  at  Monterey, 
California  is  described  and  the  oceanic  and  atmospheric  processes  which  cause 
these  variations  are  examined.  Monthly  mean  sea  level  anomalies  at  Monterey  are 
significantly  correlated  with  those  observed  at  tide  stations  as  distant  as  Prince 
Rupert,  Canada  and  Callao,  Peru,  indicating  that  these  anomalies  are  related  to 
large  scale  rather  than  to  strictly  local  atmospheric  or  oceanic  changes.  Multiple 
regression  analysis  indicates  that  monthly  anomalies  of  atmospheric  pressure,  sea 
surface  temperature,  and  meridional  wind  stress  account  for  much  of  the  observed 
monthly  sea  level  variability.  Sea  level  elevations  at  Monterey,  when  adjusted  for 
the  atmospheric  pressure  (barometric)  effect,  show  a  seasonal  high  in  winter  and  a 
minimum  in  spring.  These  seasonal  variations  are  in  phase  with  those  of  nearby 
steric  height  observations.  There  is  also  good  agreement  between  weekly  mean  sea 
level  and  steric  height  observations  in  a  time-series  sense.  Because  of  the  close 
agreement  between  sea  level  and  steric  height,  and  the  high  correlation  of 
Monterey  sea  level  with  that  at  nearby  tide  stations,  steric  height  and  sea  level 
variations  both  must  be  related  to  variations  in  the  geostrophic  current  flow. 
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I.  INTRODUCTION 


The  study  of  sea  level  and  its  changes  has  interested  man  for  many  years. 
Historical  sea  level  time-series  data  are  unique  among  marine  data  sources  in  that 
they  have  been  obtained  continuously  and  inexpensively  over  periods  of  decades  or 
longer  at  a  large  number  of  fixed  coastal  locations  worldwide.  Sea  level  records 
include  not  only  periodic  fluctuations  due  to  astronomic  tides  but  also  non-tidal  sea 
level  fluctuations  resulting  from  various  oceanic  and  atmospheric  processes.  The 
latter  can  be  isolated  as  anomalies  by  filtering  out  the  astronomic  tides,  thus 
making  measurements  of  sea  level  useful  as  a  specially  integrated  index  of 
nearshore  and  offshore  ocean  and  atmospheric  changes. 

This  paper  examines  the  character  of  sea  level  anomalies  at  Monterey, 
California  and  inquires  into  the  relative  importance  of  the  large-scale  atmospheric 
and  ocean  processes  which  affect  non-tidal  sea  level  measurements  there.  Sea 
level  variability  on  seasonal,  monthly,  weekly,  and  hourly  time  scales  is  described, 
as  are  the  physical  processes  which  cause  these  changes.  An  understanding  of 
these  processes  will  allow  the  use  of  the  long  series  of  sea  level  data  to  reconstruct 
historical  changes  in  the  oceanographic  environment  of  the  California  Current 
system,  which,  in  turn,  will  aid  in  understanding  past  changes  in  distribution, 
abundance,  and  availability  of  marine  fish  populations. 

This  study  was  supported  by  the  National  Marine  Fisheries  Service  (NMFS) 
which  is  interested  in  the  use  of  sea  level  data  for  identification  of  anomalous 
environmental  periods  and  monitoring  of  changes  in  oceanographic  conditions 
offshore. 
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A.  EARLIER  STUDIES  ON  SEA  LEVEL  VARIATIONS 

Sea  level  variations  along  the  Pacific  coast  and  their  relationship  to  various 
environmental  phenomena  have  been  examined  from  a  number  of  different  points 
of  view.  In  addition  to  the  well-understood  astronomically  induced  periodicities,  it 
is  widely  recognized  that  coastal  sea  level  measurements  may  be  influenced  by: 

1)  changes  in  atmospheric  pressure  over  the  ocean  surface, 

2)  changes  in  average  density  of  the  sea  water  column, 

3)  redistribution  of  water  mass  due  to  wind  stress, 

4)  wind  set-up  or  set-down  against  the  coast  due  to  storms, 

5)  subsidence  or  uplift  of  the  land  upon  which  the  tide  gage  is  located, 

6)  long  period  astronomic  tides, 

7)  changes  in  total  mass  of  water  in  the  oceans  associated  with  the  glacial 

ice  budget,  and 

8)  wind  waves  and  swell. 

These  physical  processes  are  discussed  by  Montgomery  (1938).  LaFond  (1939)  found 
close  agreement  between  weekly  mean  sea  level  measured  at  La  Jolla,  California 
and  the  geopotential  topography  offshore,  thus  directly  relating  ocean  currents  to 
sea  level.  Jacobs  (1939)  suggested  that  the  relationships  observed  by  LaFond  were 
not  entirely  due  to  changes  in  the  density  of  surface  water  but  rather  to  actual 
slopes  induced  by  wind-driven  water  transport  along  the  coast.  Pattullo,  et.  al. 
(1955)  found  that  south  of  40°  N  in  the  North  Pacific  Ocean,  the  seasonal  variation 
of  steric  elevation  and  sea  level  are  in  phase,  both  having  a  maximum  elevation  in 
late  summer  or  early  fall  and  a  minimum  elevation  in  winter.  This  they  took  as  a 
consequence  of  seasonal  heating  and  cooling.  These  investigators  further  found 
that  seasonal  variations  in  sea  level  north  of  40°  N  along  the  northwest  coast  of  the 
United  States  could  not  be  explained  by  steric  considerations  alone,  suggesting  that 
non-isostatic  processes  such  as  wind  and  currents  can  lead  to  appreciable  regional 
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deviations.  Roden  (1960)  used  autocorrelation  and  spectral  techniques  to  examine 
the  relationship  between  monthly  mean  sea  level  pressure,  wind,  and  sea  surface 
temperature  (SST)  at  several  stations  along  the  Pacific  coast.  He  found  good 
coherence  between  anomalies  of  sea  level  and  atmospheric  pressure,  moderate  to 
poor  coherence  between  SST  and  sea  level  depending  on  the  location  of  the  station, 
and  some  coherence  between  anomalies  of  sea  level  and  the  north-south  component 
of  the  geostrophic  wind.  Sturges  (1974)  found  high  correlations  between  occasional 
steric  observations  and  3-day  mean  sea  levels  at  Neah  Bay,  Washington  and  San 
Diego,  California.  Reid  and  Mantyla  (1976)  demonstrated  that  the  winter  increase 
in  seasonal  sea  level  elevation  along  the  northern  north  Pacific  coast  results  from 
increased  overall  flow  in  the  subarctic  cyclonic  gyre. 

B.  OCEAN  AND  ATMOSPHERIC  PROCESSES  NEAR  MONTEREY 

The  California  Current  is  a  broad,  diffuse,  southward  flowing  eastern 
boundary  current.  The  strength  of  the  current  is  affected  by  the  winds  over  the 
current  which,  in  turn,  are  controlled  by  the  strength  and  location  of  the  Aleutian 
low-pressure  cell  located  over  the  Aleutian  Islands,  the  Pacific  high  pressure  cell 
located  mainly  over  the  ocean  east  of  the  Hawaiian  Islands,  and  the  thermal  low- 
pressure  cell  located  over  the  western  United  States.  During  the  spring  and 
summer  months  the  Aleutian  low  weakens  and  the  Pacific  high  intensifies  and 
moves  northward.  Winds  over  the  current  during  this  period  are  mainly  from  the 
northwest  and  are  strongest  when  the  Pacific  high  and  thermal  low  pressure  cells 
are  closest  together  and  relatively  intense.  Winds  weaken  or  change  direction 
when  this  pressure  gradient  decreases.  The  seasonal  change  in  strength  and 
location  of  these  pressure  cells  thus  causes  seasonal  changes  in  the  winds  (Reid, 
Roden,  and  Wyllie,  1958). 
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Skogsberg  (1936)  described  three  distinct  seasonal  phases  in  his  study  of  the 
hydrography  of  Monterey  Bay.  The  calendar  year  opens  in  the  countercurrent  or 
Davidson  Current  phase.  In  late  fall  and  early  winter  of  most  years  northerly  winds 
are  weak  and  variable,  and  a  northward  flowing  countercurrent  usually  forms  at  the 
surface  along  the  central  California  coast.  This  current  is  reinforced  by  intermit¬ 
tent  periods  of  southerly  winds  during  winter.  The  general  north  northwest-south 
southeast  trend  of  the  coastline  and  the  movement  of  surface  water  to  the  right  of 
the  wind  due  to  the  Ekman  effect  cause  onshore  transport  of  surface  waters  and 
piling  up  against  the  coast.  Minimal  solar  radiation  and  strong  vertical  mixing  of 
surface  waters  by  winter  storms  decrease  sea  surface  temperatures  to  a  seasonal 
minimum  during  January  or  February.  While  SST’s  decline  during  the  Davidson 
Current  period,  temperatures  at  deeper  levels  slowly  increase  due  to  advection  of 
warm  waters  from  the  south.  For  example,  temperatures  at  50  meters  depth  reach 
a  seasonal  maximum  during  December  and  January  (Skogsberg,  1936;  Bolin  and 
Abbott,  1963).  The  end  of  the  Davidson  Current  period  is  variable  and  difficult  to 
pinpoint.  About  March,  the  offshore  high  pressure  cell  intensifies  and  northwest 
winds  become  frequent.  The  Ekman  deflection  causes  offshore  transport  of  surface 
water.  In  the  nearshore  region,  some  of  this  water  is  replaced  by  cold,  nutrient- 
rich  subsurface  water  upwelled  from  the  upper  hundred  or  so  meters.  Upwelling  is 
strongest  when  northerly  winds  are  strongest,  and  usually  reaches  a  maximum  near 
Monterey  in  May  or  June  (Bakun,  1975).  By  August,  northerly  winds  begin  to 
slacken  and  the  increased  solar  radiation  of  late  spring  and  summer  results  in  a 
steady  rise  in  SST  that  usually  continues  through  September.  September  and 
October  bring  about  a  period  of  calmer  winds  that  Skogsberg  (1936)  called  the 
oceanic  period.  With  a  slackening  of  windstress,  the  cool,  upwelled  water  begins  to 
sink  and  is  replaced  by  warmer  surface  water  from  offshore.  Coastal  SST's  rise  to 
their  highest  seasonal  values  and  strong  vertical  temperature  gradients  form  (Bolin 
and  Abbott,  1963). 
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To  summarize,  the  oceanographic  regime  off  Monterey  is  marked  by  three 
distinct  periods:  The  Davidson  Current  period,  occurring  during  November  through 
February,  has  weak  northerly  winds,  strong  winter  storm  events,  northward  current 
flow,  and  onshore  transport  of  surface  water.  The  upwelling  period,  occurring  in 
March  through  August,  has  strong  northwest  winds,  southward  current  flow, 
offshore  transport  of  surface  water,  and  upwelling  of  cool,  nutrient-rich  water. 
The  oceanic  period,  occurring  during  September  and  October,  is  a  period  of  calm 
between  the  northerly  winds  of  the  upwelling  period  and  the  southerly  winds  of 
winter.  During  this  period,  surface  temperatures  increase  and  strong  vertical 
temperature  gradients  form.  These  are  the  average  seasonal  characteristics  in  the 
meteorological  and  oceanic  regimes  affecting  Monterey,  however,  there  are 
marked  year-to-year  differences  in  both  timing  and  intensity  of  the  events 
described. 


n.  DESCRIPTION  OF  DATA 

Recorded  tide  data  from  the  tide  station  at  Monterey,  California  were  chosen 
for  analysis  because  the  gage  lies  in  a  coastal  area  of  interest  to  the  NMFS  and  is 
exposed  to  open  ocean  conditions  with  no  nearby  river  discharge  to  affect  sea  level 
measurements.  The  Monterey  gage  is  the  only  primary  tide  station  maintained  by 
the  National  Ocean  Survey  (NOS)  between  San  Francisco  and  Avila,  and  thus  fills  a 
large  data  gap  along  the  central  California  coast.  The  Monterey  station  has  been 
operated  continuously  since  1983  by  the  Naval  Postgraduate  School  under  the 
direction  of  Dr.  Warren  C.  Thompson  but  the  time-series  data  were  not  previously 
fully  analyzed. 

Monterey  Bay  is  located  about  120  km  south  of  San  Francisco,  California. 
The  bay,  which  is  bisected  by  a  deep  submarine  canyon,  is  a  large,  semi-elliptical 
coastal  feature  measuring  about  37  km  wide  at  the  mouth  and  about  19  km  from  the 
mouth  to  the  innermost  point.  The  tide  station  is  located  along  the  southern  edge 
of  the  bay  near  the  end  of  Monterey  Municipal  Wharf  #2  in  a  water  depth  of 
approximately  6.8  m.  Because  of  the  open  shape  of  the  bay  and  the  narrow  width 
of  the  continental  shelf,  tide  measurements  obtained  here  are  presumed  to  fully 
represent  those  of  the  open  coast. 

In  addition  to  sea  level  data,  meterological  and  oceanographic  data  repre¬ 
sentative  of  the  Monterey  area,  including  surface  atmospheric  pressure  data, 
geostrophic  wind  data,  surface  salinity  and  temperature  data,  and  deep  hydrocast 
data  were  used  in  this  study.  The  general  proximity  of  the  various  data  sources 
allowed  direct  comparison  of  variables  with  minimal  problems  resulting  from 
spatial  distortion.  Figure  1  shows  the  location  from  which  each  of  the  data  sources 
were  derived,  along  with  the  nearby  bathymetry. 
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A.  MONTEREY  SEA  LEVEL  DATA 

1.  Tide  Gages 

A  standard  recording  tide  gage,  which  traces  tide  heights  continuously  on  a 
strip  chart,  was  installed  at  the  Monterey  tide  station  by  NPS  personnel  in  June 
1963.  This  analog  system  is  entirely  mechanical  and  is  highly  dependable  when 
maintained  properly.  A  drum-mounted  strip  chart  is  rotated  by  a  spring-driven 
clock  mechanism,  and  a  pencil  records  sea  level  changes  by  means  of  a  float-pulley 
system.  A  second  instrument,  a  Fischer-Porter  digital  tide  gage,  was  installed 
adjacent  to  the  analog  gage  by  the  National  Ocean  Survey  in  November  1973.  This 
is  an  electrically  operated  system  which  punches  digital  data  on  foil  tape.  Both 
gages  use  the  same  21.6  cm  diameter  float  and  have  operated  simultaneously  since 
November  1973.  The  stilling  well,  which  serves  as  a  low  pass  filter  for  oscillations 
with  periods  greater  than  a  minute,  consists  of  a  30.5  cm  diameter  steel  pipe  with 
a  2.5  cm  diameter  orifice  at  the  bottom.  Both  gages  are  checked  for  accuracy  of 
time  and  height  and  are  annotated  about  five  times  per  week. 

2,  Data  Processing  and  Reduction 

Continuous  tide  traces  obtained  from  the  analog  gage  during  the  period  July 
20,  1963  through  December  31,  1974  were  manually  digitized  for  use  in  this  study. 
Digitization  was  performed  by  Ocean  Data  Systems,  Inc.,  Monterey,  CA  under 
contract  to  the  NMFS  and  with  NPS  guidance.  Datums  were  reviewed  and  data 
were  reduced  to  hourly  sea  level  heights  using  standard  NOS  procedures  (Coast  and 
Geodetic  Survey,  1965).  Data  from  the  digital  gage  for  the  period  January  1,  1974 
through  September  31,  1976  were  processed  for  hourly  heights  by  the  NOS  and 
provided  for  use  in  this  study.  Data  from  both  gages  were  recorded  in  feet  and 
later  converted  to  centimeters  in  this  study.  The  hourly  heights  are  resolved  to  the 
nearest  0.1  foot  (3.0  cm)  and  times  of  observation  (Pacific  Standard  Time)  are 
accurate  to  within  six  minutes.  A  small  percentage  of  the  hourly  sea  level  data 

were  missing,  either  rejected  as  erroneous  or  lost  due  to  equipment  malfunctions. 
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As  a  result  some  monthly  means  contain  less  than  a  full  month  of  data.  Missing 
data  of  duration  of  a  day  or  longer  are  listed  in  Appendix  A. 


All  hourly  heights  were  measured  relative  to  the  station  datum  established  by 
the  NOS  in  November,  1973.  Mean  sea  level  for  the  period  1963  through  1978  lies  at 
184.4  cm  and  the  National  Geodetic  Vertical  Datum  lies  182.88  cm  above  the 
station  datum. 

3.  Merging  of  Analog  and  Digital  Tide  Data 

To  obtain  the  longest  possible  continuous  tide  record,  it  was  necessary  to 
merge  the  analog  and  the  more  recent  digital  data  sets.  Before  the  data  sets  were 
combined,  the  response  of  the  two  gages  was  analyzed  by  comparing  the  hourly 
heights  from  both  tide  records  for  the  calendar  year  1974.  The  correlation 
coefficient  between  the  analog  and  digital  data  sets  exceeds  0.99,  as  anticipated. 
The  regression  equation  for  the  two  sets  of  hourly  heights  is  Y  =  1.48  cm  +  0.98X, 
where  X  refers  to  digital  data  and  Y  refers  to  analog  data  in  cm. 

Figure  2  shows  a  histogram  of  the  difference  (digital  minus  analog  data) 
between  the  two  sets  of  hourly  sea  levels  for  the  calendar  year  1974.  The  mean 
difference  was  found  to  be  -0.06  cm.  The  frequency  distribution  of  the  differences 
resembles  a  normal  distribution,  with  a  standard  deviation  of  3.7  cm.  Nearly  all  of 
the  differences  are  attributed  to  the  fact  that  the  digital  data  were  recorded  as 
instantaneous  values,  which  can  include  short-term  sea  level  fluctuations  such  as 
long  waves  or  seiches,  whereas  in  the  analog  data  these  short-term  fluctuations 
were  filtered  out  by  manually  smoothing  the  tide  curve  before  digitizing. 

It  was  concluded  that  differences  between  the  two  data  sets  were  negligible, 
and  that  the  analog  and  digital  data  could  be  combined  without  bias.  Thus,  analog 
data  from  the  period  July  20,  1963  through  December  31,  1974  were  combined  with 
digital  data  from  the  period  January  1,  1975  through  August  31,  1976  to  form  a  13- 
year  time  series  containing  a  total  of  107,954  hourly  observations. 


Relative  Frequency  (%) 


4.  Long  Period  Sea  Level  Changes 


Tide  gages  show  the  height  of  the  sea  level  relative  to  land.  Thus,  changes  in 
mean  sea  level  over  periods  of  years  or  decades  can  result  from  the  addition  or 
removal  of  water  from  the  oceans  due  to  global  climatic  variations,  from 
subsidence  or  emergence  of  the  land  upon  which  the  gage  is  located,  or  from  long- 
period  astronomic  tides.  For  example,  some  long-period  trends  in  sea  level 
records,  such  as  the  rise  in  sea  level  in  Panama  described  by  Roden  (1963)  or  the 
drop  in  sea  level  in  the  Juneau,  Alaska  area  described  by  Hicks  (1973)  clearly  result 
from  local  or  regional  land  subsidence  or  uplift. 

To  determine  trends  in  the  Monterey  sea  level  record  during  the  period  1963 
through  1978,  the  time-series  of  monthly  means  was  analyzed  using  a  least-squares 
fit.  This  analysis  yielded  a  relative  rise  in  sea  level  of  0.01  cm/year.  The 
variability  in  sea  level  due  to  oceanographic  and  meterological  processes  thus 
greatly  exceeds  measured  trends.  Accordingly,  the  effects  of  trends  were 
neglected  in  this  study. 

Of  the  long  period  astronomic  tides,  the  nodal  tidal  constituent,  which  results 
from  the  changing  declination  of  the  moon  over  a  period  of  18.61  years,  has  the 
greatest  amplitude.  The  theoretical  amplitude  of  this  constituent  varies  with 
latitude,  with  maximum  effects  at  the  equator  and  the  poles  and  minimum  effects 
near  the  latitudes  35°  N  and  35°  S  (Lisitzin,  1974).  A  second  significant  long  period 
constituent,  the  annual  solar  tide,  has  an  amplitude  approximately  one  fifth  of  the 
nodal  tide  component.  The  effects  of  this  tidal  constituent  vary  with  latitude  in  a 
manner  similar  to  that  of  the  nodal  tide.  Monterey,  located  near  36°  N,  is  in  a 
region  where  the  ranges  of  both  of  these  long  period  tides  are  about  1  cm,  so  these 
effects  were  neglected  in  this  study. 


B.  OCEAN  AND  ATMOSPHERIC  DATA 

The  computed  atmospheric  pressure  and  wind  data  used  in  this  study  were 
derived  from  six-hourly  synoptic  surface  pressure  fields  prepared  by  Fleet 
Numerical  Oceanography  Center  (FNOC).  The  pressure  fields,  arranged  on  a  grid 
with  a  mesh  length  of  3  degrees  latitude,  were  used  to  compute  geostrophic  winds, 
from  which  windstress,  Ekman  transport,  and  Sverdrup  transport  estimates  were 
calculated  at  a  deep  water  site  approximately  14  km  west  of  Monterey  (Figure  1). 
A  complete  description  of  the  methods  and  computations  used  in  these  calculations 
is  given  by  Bakun  (1975).  Briefly,  the  geostrophic  wind  was  computed  at  the  point 
36.6°  N,  122.1°  W  and  an  estimate  of  the  wind  near  the  sea  surface  was  made  by 
rotating  the  geostrophic  wind  vector  15  degrees  to  the  left  and  reducing  its 
magnitude  by  30%.  The  surface  wind  stress  was  computed  and  the  wind  stress 
vector  was  resolved  into  north-south  (meridional  or  alongcoast)  and  east-west 
(zonal  or  acrosscoast)  components.  Ekman  transport  was  computed  and  offshore- 
onshore  transport  was  determined  by  resolving  the  vector  component  perpendicular 
to  the  general  trend  of  the  coastline.  Sverdrup  transport  was  calculated  as 
described  by  Nelson  (1977). 

The  surface  temperature  and  salinity  data  were  obtained  from  samples  taken 
daily  at  Hopkins  Marine  Station  of  Stanford  University  during  the  period  January 
1963  to  May  1975  (Scripps  Institution  of  Oceanography,  1963  to  1976).  SST  data  from 
June  1975  to  December  1978  were  taken  at  the  Monterey  tide  station  by  Mr.  Jerry 
Norton  of  the  Naval  Postgraduate  School.  Salinity  data  from  Hopkins  are  not 
available  later  than  May  1975. 

Monthly  means  and  anomalies  of  sea  level,  and  of  the  ocean  and  atmospheric 
data  described  in  the  above  sections  are  presented  graphically  and  in  tabular  form 
in  Appendix  B. 

To  examine  the  relationship  between  sea  level  and  dynamic  height,  a  series  of 
hydrographic  cast  data  were  assembled  for  a  station  located  in  mid-Monterey  Bay, 
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about  19  km  northwest  of  the  tide  station  (Figure  1).  This  hydrographic  station  is 
located  near  the  mouth  of  the  Monterey  submarine  canyon  where  the  water  depth 
is  over  900  meters.  The  hydrographic  cast  data  were  taken  semi-monthly  by 
Hopkins  Marine  Station.  Sampling  during  the  first  years  of  the  program  was 
limited  to  the  upper  50  meters  of  the  water  column  but  in  1968  the  sampling  depth 
was  increased  to  over  500  meters  (Hopkins  Marine  Station,  1968-1973).  Sampling 
was  discontinued  by  Hopkins  in  December  1973  and  was  resumed  on  a  semi-monthly 
basis  by  Moss  Landing  Marine  Laboratory  from  July  1974  to  June  1978  (Broenkow,  et 
al.,  1975  and  1976;  Lesley,  1977;  Chinburg,  et  al.,  1978).  Hydrographic  data  for  the 
ten-year  period  January  1968  through  December  1977  were  digitized,  and  long-term 
monthly  mean  dynamic  heights  and  dynamic  height  anomalies  were  calculated. 
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III.  SEA  LEVEL  AT  MONTEREY 


The  tides  at  Monterey  are  mixed,  predominantly  semidiurnal,  and  are 
composed  of  two  high  and  two  low  waters  per  24.8  hour  tidal  cycle  (Figure  3). 
Analysis  of  the  13  years  of  hourly  data  show  that  the  largest  daily  tide  range 
recorded,  274.3  cm,  occurred  on  December  20,  1968.  The  maximum  water  level 
during  the  period  of  record,  329.2  cm  above  station  datum,  occurred  on  January  18, 
1973  and  was  30.5  cm  above  the  predicted  tide  for  that  time.  The  minimum  water 
level,  36.6  cm  above  station  datum,  occurred  on  June  11,  1968  and  was  15.2  cm 
below  the  predicted  tide.  This  analysis  does  not  include  data  for  March  28,  1964 
because  of  a  tsunami  resulting  from  the  Alaskan  earthquake. 

It  is  recognized  that  the  time  series  of  hourly  sea  levels  could  be  analyzed  for 
the  occurrence,  amplitude,  and  duration  of  anomalous  sea  level  events.  This  was 
done  by  Maixner  (1973)  who  examined  Monterey  sea  level  anomalies  for  the  year 
1971.  It  was  decided  for  this  study,  however,  to  concentrate  on  variations  identified 
in  mean  monthly  sea  level  data  and  on  their  atmospheric  and  oceanographic  causes; 
weekly  mean  and  six-hourly  sea  level  data  are  also  examined  in  a  limited  way.  The 
statistical  characteristics  of  hourly  deviations  from  the  predicted  sea  level  are  also 
examined. 

A.  MEANS  AND  VARIATIONS 

1.  Hourly  Sea  Level 

To  analyze  non-tidal  sea  level  variations,  which  are  small  compared  to  the 
normal  tide  range  in  this  area,  the  tidal  signal  must  be  removed  from  the  data. 
This  can  be  done  by  averaging,  filtering,  or  subtracting  predicted  tides  from  the 
data.  The  latter  method  was  used  in  this  study.  The  Tide  Predictions  Branch  of 
the  NOS  performed  a  harmonic  analysis  of  365  days  of  hourly  Monterey  tide  heights 
and  isolated  37  harmonic  constituents.  These  constituents  are  listed  by  Maixner 


(1973).  Using  the  20  constituents  whose  amplitudes  were  greater  than  0.61  cm,  the 
NOS  computed  predicted  hourly  tide  heights  for  the  period  of  record.  Predicted 
hourly  heights  were  then  subtracted  from  the  13  years  of  observed  hourly  heights  to 
yield  non-astronomic  residuals.  The  frequency  of  occurrence  of  these  sea  level 
differences  (observed  minus  predicted),  which  total  nearly  108,000  values,  approxi¬ 
mates  a  normal  or  Gaussian  distribution  (Figure  4).  94.5%  of  the  observations  lie 
within  15.2  cm  (0.5  foot)  of  the  predicted  tide  and  99.9%  lie  within  30.5  cm  (1.0 
foot).  The  maximum  observed  difference  was  39.6  cm.  The  standard  deviation  of 
the  differences  was  8.7  cm,  skewness  -0.02,  and  kurtosis  3.2. 

The  distribution  of  hourly  differences  describes  non-tidal  sea  level  variations 
over  a  13  year  period  but  gives  no  information  about  seasonal  variation.  Does  the 
frequency  distribution  change  from  month  to  month?  Are  distributions  for  winter 
months  the  same  as  those  for  summer?  To  define  the  seasonal  change,  curves  were 
generated  using  data  from  8,200  to  9,800  observations  for  each  of  the  12  months  of 
the  year  and  these  are  shown  in  Figure  5.  The  frequency  distribution  of  non-tidal 
sea  level  fluctuations  is  seen  in  the  figure  to  change  seasonally.  In  April,  for 
example,  73%  of  the  observed  sea  levels  were  lower  than  predicted,  but  in 
September,  81%  of  the  observed  data  were  greater  than  predicted.  From  March 
thru  May,  observed  sea  levels  tend  to  be  lower  than  predicted  sea  levels,  probably 
due  to  offshore  Ekman  transport,  SST,  and  atmospheric  pressure  effects  as 
discussed  in  another  section.  From  July  through  January,  observed  sea  levels  are 
greater  than  predicted  due  to  atmospheric  pressure  and  thermal  effects  during 
summer  and  fall,  and  to  onshore  transport,  pressure,  and  thermal  effects  during  the 
Davidson  Current  Period  in  December  and  January. 

The  distributions  for  winter  months  are  wider  and  less  peaked  than  those  of 
summer  months,  indicating  greater  variability  and  larger  non-tidal  events.  The 
distributions  for  July  and  August  are  notably  narrow  and  peaked  in  contrast. 
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•igure  4.  Difference  Between  Observed  and  Predicted  Hourly  Tide 
Heights  (1963-1976) 
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2.  Monthly  Mean  Sea  Level 

Averaging  over  intervals  of  weeks  to  months  removes  the  effects  of  the 
principal  diurnal,  semi-diurnal,  and  other  short-term  tidal  components  from  the 
data,  reduces  large  quantities  of  data  to  manageable  size,  and  emphasizes  time 
scales  important  for  many  marine  biological  processes. 

The  monthly  means  were  calculated  for  the  period  July  1963  through  August 
1976  by  the  author,  and  for  the  period  September  1976  through  December  1978  by 
the  NOS.  Figure  6  shows  the  long  term  monthly  means,  standard  deviations,  and 

extreme  monthly  sea  levels  based  on  the  16-year  period  1963  to  1978.  Mean  sea 

level  is  seen  to  be  lowest  in  April  and  highest  in  September,  with  a  mean  annual 
range  of  13.6  cm.  Variability  is  highest  during  winter  months,  with  monthly 
standard  deviations  during  winter  being  almost  double  those  for  summer.  The 
range  between  maximum  and  minimum  monthly  values  reaches  a  high  of  21.0  cm  in 
January  and  a  low  of  8.5  cm  in  August. 

Monthly  sea  level  anomalies  were  calculated  as  differences  between  the 

monthly  mean  and  the  long-term  mean  for  the  same  month.  Calculation  of 

anomalies  in  this  manner  removes  the  annual  cycle  from  the  data  and  allows  direct 
comparison  of  month-to-month  variability.  Monthly  mean  sea  levels  and  their 
anomalies  are  shown  in  tabular  and  graphical  form  in  Figures  7A-7C.  In  these 
figures  extreme  monthly  sea  level  anomalies  are  shown  to  range  from  -10.8  cm  in 
December  1975  to  +10.7  cm  in  January  1978.  Periods  of  anomalously  high  sea  level 
occurred  during  1969,  1972-1973,  1976-1977,  and  early  1978,  and  periods  of  anoma¬ 
lously  low  sea  level  occurred  in  1964,  1970,  1971,  1973,  1975-1976,  and  1977. 

To  statistically  define  the  persistence  of  anomalous  periods,  the  auto¬ 
correlation  function  was  used.  This  function  describes  the  decay  of  the  correlation 
coefficient  as  the  data  series  is  time  shifted  relative  to  itself  an  increasing  number 
of  months.  The  autocorrelation  function  of  monthly  Monterey  sea  level  anomalies, 
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Figure  7A.  Monterey  Monthly  Sea  Level  Means  and  Anomalies 


(Monthly  Means  are  Shown  as  Heavy  Lines  and  Anomalies 
as  Light  Lines) 
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Figure  7B. 
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Figure  7C. 


Figure  8,  shows  sea  level  anomalies  to  be  correlated  at  the  5%  level  of  significance 
for  lags  of  up  to  5  months,  indicating  that  anomalies  persist  over  a  period  of 
several  months.  The  autocorrelation  function  of  the  sea  level  series  appears  to 
decay  exponentially  for  the  first  8  months  or  so,  with  significant  negative 
autocorrelation  coefficients  occuring  from  lags  of  11  to  18  and  23  to  26  months. 

B.  RELATION  TO  OTHER  PACIFIC  COAST  TIDE  STATIONS 

We  have  seen  that  mean  monthly  sea  level  anomalies  at  Monterey  tend  to 
persist  for  about  5  months.  The  question  naturally  arises  as  to  whether  these 
anomalies  are  of  local  or  regional  geographic  extent.  To  determine  the  spacial  and 
temporal  coherence  between  the  monthly  anomalies  at  Monterey  and  those 
observed  at  neighboring  tide  recording  stations,  monthly  mean  data  were  assembled 
for  15  tide  stations  along  the  Pacific  coast  ranging  from  Sitka,  Alaska  to  Callao, 
Peru  (Figure  9).  These  data  were  obtained  from  Dr.  Klaus  Wyrtki  of  the  University 
of  Hawaii  and  from  the  NOS.  Stations  selected  for  analysis  were  those  having  the 
best  combination  of  the  following  characteristics:  Representativeness  of  open 
ocean  conditions,  a  long  and  continuous  data  record,  a  constant  tidal  reference 
datum  during  the  time  period  of  interest,  and  suitable  spacing  between  station 
locations  along  the  coast.  For  each  station  long-term  monthly  means  were 
calculated  from  the  available  tide  measurements  for  the  period  1963  to  1978  and 
monthly  sea  level  anomalies  were  derived,  as  for  Monterey.  These  data  are  shown 
in  Figures  10A-10C. 

Variation  in  the  month-to-month  value  of  the  anomaly  may  be  seen  from  the 
figures  to  be  greatest  for  stations  north  of  Crescent  City,  showing  the  effects  of 
energetic  winter  storms.  Perhaps  the  most  striking  feature  of  the  time  series  is 
the  high  visual  correlation  of  anomalies  along  the  coast.  The  periods  of 
anomalously  high  sea  level  at  Monterey  during  1969,  1972-1973,  1976-1977,  and  1978 
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DATA  POINTS  180 

NUMBER  OF  LAGS  40 

MEAN  -.262 

ST  DEVIATION  -3.996 

VARIANCE  15.972 

SKEWNESS  -.030 

KURTOSIS  3.170 


Figure  8.  Autocorrelation  Function  for  Monthly  Sea  Level  at  Monterey 
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were  common  to  most  stations  where  data  are  available.  Similarly,  the  periods  of 
anomalously  low  sea  level  seen  at  Monterey  in  1964,  1970,  1971,  1973,  1975-1976,  and 
1977  occurred  at  most  of  the  other  stations.  Also,  the  anomalies  tend  to  reverse 
their  sign  more  frequently  at  the  northern  stations  than  at  the  southern  stations. 
This  observation  appears  worthy  of  future  study. 

Correlations  of  the  monthly  sea  level  anomalies  between  stations  were 
calculated  using  the  BMDP8D  statistical  program  (Dixon,  1975)  and  are  tabulated  in 
Table  1.  The  correlation  of  the  selected  tide  stations  relative  to  Monterey  is  shown 
graphically  in  Figure  11.  Correlation  of  the  Monterey  anomalies  is  seen  to  be 
highest  with  San  Francisco  (r  =  0.85)  and  lowest  with  Sitka  (r  =  0.15).  It  can  also  be 
seen  that  the  correlation  coefficient  drops  off  more  rapidly  with  distance  to  the 
north  of  Monterey  them  to  the  south. 

Osmer  and  Huyer  (1978)  suggested  the  existence  of  two  domains  of  coastal 
sea  level  fluctuations,  with  a  boundary  located  south  of  San  Francisco  in  winter  and 
north  of  Crescent  City  in  the  spring  and  summer.  The  general  location  of  their 
break-point  is  in  agreement  with  the  findings  of  Zee  (1975),  who  suggested  that  sea 
level  anomalies  in  his  southern  group  of  stations  were  related  to  non-seasonal 
vertical  movement  of  the  thermocline.  That  a  strong  gradient  or  boundary  may 
exist  between  northern  and  southern  stations  is  further  suggested  by  Nelson  (1977) 
who  showed  that  the  area  off  northern  California  near  Cape  Mendocino  is  one  of 
marked  change  in  the  seasonal  surface  wind  stress  field.  The  mean  seasonal  wind 
stress  field  over  the  coastal  ocean  south  of  Cape  Mendocino  is  alongshore 
(southward)  all  year  while  the  stress  field  north  of  Cape  Mendocino  is  strongly 
onshore  in  winter  and  alongshore  (southward)  in  summer. 

The  alongcoast  extent  of  sea  level  anomalies  observed  at  Monterey  was 
further  examined  in  a  time-distance  domain.  The  monthly  anomalies  from  the 
series  of  15  coastal  stations  from  Sitka,  Alaska  to  Callao,  Peru  were  plotted  and 
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11.  Correlation  of  Monthly  Sea  Level  Anomalies  at  Selected 
West  Coast  Tide  Stations  Relative  to  Monterey 
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contoured  at  5  cm  intervals  for  the  period  1963  to  1974  (Figures  12A-12D).  Data  for 
the  years  1975  to  1978  were  not  available  for  several  stations  so  plots  for  these 
years  are  not  included.  Anomalies  are  seen  to  fall  into  recognizable  patterns  which 
are  coherent  in  both  time  and  space.  In  Figure  12A,  for  example,  large  negative 
anomalies  can  be  seen  in  January  1963  extending  from  Crescent  City  to  Sitka  and 
large  positive  anomalies  in  the  same  region  occur  in  the  subsequent  fall  and  winter. 

Larger  anomalies  and  larger  anomaly  gradients  occur  northward  of  a 
boundary  lying  generally  between  Crescent  City  and  Monterey.  Anomalous  events 
to  the  north  of  this  boundary  tend  to  occur  simultaneously  along  the  coast  and  are 
persistent  for  one  or  two  months.  Anomaly  magnitudes  and  gradients  are  also 
generally  larger  southward  of  a  second  diffuse  boundary  zone  lying  approximately 
between  Manzanillo  and  Quepos.  Between  these  zones  the  anomaly  field  is 
relatively  flat.  Southward  of  the  general  boundary  between  Crescent  City  and 
Monterey  sea  level  anomalies  are  of  longer  duration,  as  was  noted  earlier  in 
reference  to  Figure  10. 

A  particularly  interesting  event  is  the  anomalously  high  sea  level  during  the 
period  October  1972  through  February  1973  between  Callao  and  San  Francisco 
(Figure  12D).  This  was  a  period  of  strong  El  Nino  activity  in  the  eastern  tropical 
Pacific.  During  El  NilTo  occurrences  there  is  a  rapid  rise  in  sea  level  in  the  eastern 
tropical  Pacific  accompanied  by  a  fall  in  sea  level  in  the  western  Pacific  (Wyrtki, 
1977).  In  Figure  12D,  for  example,  a  peak  anomaly  of  25  cm  was  observed  at 
Manzanillo  in  December  1972,  where  the  occurrence  of  high  sea  levels  preceeded 
those  observed  at  more  northern  stations  by  a  month  or  more.  At  Monterey,  sea 
levels  were  higher  than  average  during  the  winter  of  1972-1973  (see  also  Figure  10). 
During  the  El  Nino  period,  as  shown  in  Appendix  B,  atmospheric  pressures  at 
Monterey  were  less  than  average  and  wind  stress  was  negligible  (except  during 
February  1973  when  anomalous  southerly  winds  resulted  in  onshore  transport  of 


(Contour  Interval  is  5  cm) 


surface  waters  and  downwelling).  Sea  surface  temperatures  were  also  anomalously 
high  from  August  1972  to  March  1973. 

In  summary,  monthly  sea  level  anomalies  at  Monterey  are  related  to  large- 
scale  influences  rather  than  to  strictly  local  events.  Table  1  shows  that  the 
anomalies  at  Monterey  are  correlated,  at  the  5%  level  of  significance,  with 
anomalies  recorded  at  stations  from  Prince  Rupert,  Canada  to  Callao,  Peru,  but 
are  more  closely  related  to  events  affecting  sea  levels  in  the  group  of  stations 
from  Crescent  City  to  Quepos.  Processes  producing  the  El  Nino  phenomenon  along 
the  coast  of  Peru  also  apparently  affect  sea  level  at  Monterey. 
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IV.  CAUSES  OF  SEA  LEVEL  VARIATIONS  AT  MONTEREY 


The  effects  of  changes  in  atmospheric  pressure,  changes  in  water  mass 
characteristics  due  to  changes  in  alongcoast  currents,  and  changes  in  average 
density  of  the  water  column  on  sea  level  are  all  interrelated.  A  change  in  the 
distribution  of  atmospheric  pressure  over  the  ocean  surface  will  generally  change 
the  horizontal  gradient  of  pressure,  resulting  in  a  change  in  the  geostrophic  and 
other  wind  components,  and  thus  in  wind  stress.  A  change  in  wind  stress  will 
change  the  wind-driven  current,  redistribute  the  mass,  and  change  the  average 
density  of  the  water  column.  W  ind  stress  changes  also  alter  wind-induced  set-up  or 
set-down  against  the  coast.  All  of  these  processes  combine  to  affect  sea  level. 

So  as  to  separate  these  effects,  correlation,  regression,  and  spectral  analysis 
techniques  were  used.  Fluctuations  in  sea  level  and  other  variables  occur  on 
various  time  scales.  For  this  reason  the  following  section  is  organized  generally  by 
time-sampling  and  specifically  by  analysis  proceedures  used.  The  procedures  used 
were  chosen  as  appropriate  for  the  character  of  the  data  to  be  analyzed. 


A.  CORRELATION  ANALYSIS 

Mean  monthly  anomalies  for  the  period  1963-1978  were  calculated  for  the 
following  oceanic  and  atmospheric  variables:  Surface  atmospheric  pressure, 
meridional  wind  stress,  zoned  wind  stress,  offshore  Ekman  transport,  Sverdrup 
transport,  SST,  and  salinity.  These  data  are  presented  in  Appendix  B.  Correlations 
between  these  variables  and  the  monthly  sea  level  anomalies  at  Monterey  were 
calculated  using  the  BMDP8D  statistical  program  and  the  results  given  in  Table  2. 
The  correlation  analysis  measures  the  strength  of  the  linear  relationships  between 
independent,  random  variables.  However,  the  variables  dealt  with  here  are  neither 
random  nor  independent  so  some  care  must  be  used  in  interpretation  of  the 
statistical  results.  In  the  following  paragraphs  each  variable  will  be  treated  in  turn 
and  the  results  of  the  correlation  analysis  discussed. 


Table  2.  Inter-Correlation  of  Monthly  Sea  Level  Anomalies 
and  Monthly  Anomalies  of  Various  Ocean  and  Atmos¬ 
pheric  Variables 
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SL 

SL 

PRES 

WS 

WS 

TSPT 

TSPT 

SAL 

SL 

1.00 

ASL 

.95 

1.00 

PRES 

-.69 

-.46 

1.00 

MWS 

.43 

.41 

-.28 

1.00 

zws 

C-.13) 

-.18 

(-.03) 

-.47 

1.00 

EKM  TSPT 

-.42 

-.41 

.25 

-.99 

.58 

1.00 

SVP  TSPT 

C.00) 

(-.06) 

-.18 

-.32 

.14 

.32 

1.00 

SAL 

-.35 

-.30 

.29 

-.31 

(.07) 

.30 

.20 

1.00 

SST 

.61 

.64 

-28 

.37 

-.17 

-.37 

(-.05) 

-.37 

Correlation  coefficients  enclosed  in  parentheses  are  not 
significant  at  the  SI  level 
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The  effect  on  sea  level  of  changes  in  atmospheric  pressure  over  the  oceans 
has  been  examined  by  a  number  of  authors  (Patullo,  et  al.,  1955;  Saur,  1962;  Roden, 
1960).  An  increase  (decrease)  in  atmospheric  pressure  results  in  a  decrease 
(increase)  in  sea  level.  These  effects  can  be  quite  large  in  some  areas,  particularly 
in  the  Gulf  of  Alaska  where  winter  storms  sire  intense  or  along  the  Gulf  or  Atlantic 
coasts  of  the  United  States  during  the  passage  of  hurricanes. 

The  isostatic  contribution  of  atmospheric  pressure  variations  to  variations  in 
sea  level  is  computed  from  the  hydrostatic  equation,  p  =  -dgh  where  p  is  the 

change  in  atmospheric  pressure  in  millibars  (mb),  d  is  the  density  of  water  in 

3  2 

g/cm  ,  g  is  the  acceleration  of  gravity  in  cm/sec  ,  and  h  is  the  change  in  sea  level 

3 

in  cm.  Applying  this  equation  to  sea  water  of  density  1.025  g/cm  and  using  980.7 
2 

cm/sec  as  the  acceleration  of  gravity,  we  find  that  an  increase  in  atmospheric 
pressure  of  1  mb  will  result  in  a  0.995  cm  depression  of  sea  level. 

The  seasonal  range  of  monthly  mean  atmospheric  pressure  at  Monterey  during 
the  period  1963  to  1978  was  7.3  mb,  but  pressure  changes  several  times  greater  than 
this  are  not  uncommon  during  the  passage  of  intense  winter  storms.  Thus,  the 
effect  of  atmospheric  pressure  is  expected  to  account  for  a  significant  portion  of 
sea  level  variability  near  Monterey. 

Maixner  (1973)  examined  hourly  data  recorded  from  the  Monterey  tide  gage 
during  the  year  1971  and  concluded  that  hourly  sea  level  responds  to  pressure 
changes  in  an  approximately  hydrostatic  manner.  The  coefficient  of  correlation 
between  monthly  mean  sea  level  anomalies  and  pressure  anomalies,  based  on  180 
months  of  simultaneous  data  from  the  period  July  1963  through  December  1978,  was 
found  in  the  present  study  to  be  -0.69  (Table  2).  The  relatively  large  negative 
correlation  coefficient  indicates  that  hydrostatic  equilibrium  is  somewhat  appli¬ 
cable  to  monthly  statistics. 
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It  was  considered  desirable  to  remove  the  static  effects  of  atmospheric 
pressure  from  the  monthly  sea  level  data  so  that  the  influence  of  other  variables  on 
sea  level  could  more  readily  be  examined.  To  accomplish  this,  monthly  mean  sea 
levels  were  adjusted  for  monthly  pressure  effects  by  increasing  (decreasing)  sea 
level  1.00  cm  for  every  1.00  mb  increase  (decrease)  of  atmospheric  pressure.  The 
use  of  the  more  accurate  value  of  0.995  was  not  warrented  in  this  study.  The 
magnitude  of  the  pressure  correction  was  determined  by  subtracting  the  long  term 
mean  pressure  for  the  period  January  1963  through  December  1978  (1016.85  mb) 
from  the  monthly  mean  atmospheric  pressures.  Mean  monthly  sea  levels  and  sea 
level  anomalies  from  which  the  hydrostatic  effect  associated  with  monthly 
pressure  anomalies  have  been  removed  are  referred  to  in  this  study  as  adjusted  sea 
levels.  The  time  series  of  adjusted  and  unadjusted  sea  level  anomalies  and  long¬ 
term  means  for  the  period  1963  through  1978  are  shown  in  Figure  13. 

In  general,  the  contribution  of  atmospheric  pressure  is  seen  in  the  time  series 
to  be  small  compared  with  the  observed  departures  of  sea  level.  In  most  months 
the  pressure  correction  is  opposite  in  sign  to  the  sea  level  anomaly  and  causes  a 
reduction  in  the  sea  level  variability.  The  greatest  differences  occur  in  winter 
months.  The  effect  of  the  subtraction  of  static  pressure  on  the  seasonal  sea  level 
shown  in  the  figure  is  to  reduce  the  range  of  the  monthly  values,  and  to  a  very 
small  extent  the  seasonal  range,  but  also  to  shift  the  time  of  highest  sea  level  from 
September  to  December.  Pressure  effects  account  for  a  portion  of  the  sea  level 
variability  but  significant  non-barometric  residuals  remain,  indicating  the  effects 
of  dynamic  as  well  as  static  processes. 

The  effects  of  wind  stress  on  sea  level  are  twofold;  the  direct  elevation  or 
depression  of  water  by  winds  normal  to  the  coast  and  the  sea  surface  slopes 
created  by  offshore  or  onshore  Ekman  transport  produced  by  winds  parallel  to  the 
coast.  The  direct  piling  up  of  water  is  commonly  observed  in  areas  with  wide, 


50 


shallow  continental  shelves  op  long,  narrow  embayments.  The  magnitude  of  this 
effect  is  dependent  on  basin  configuration,  surface  wind  velocity,  and  the  depth  of 
water.  The  shelf  in  the  Monterey  area  is  quite  narrow  with  deep  water  located 
close  inshore  so  that  the  effects  of  wind  set-up  are  small.  Defant  (1961)  showed, 
for  example,  that  a  constant  10  m/s  wind  blowing  over  a  basin  50  m  deep  would 
produce  a  sea  surface  slope  of  6.6  cm/100  km.  The  50  m  contour  near  Monterey  is 
less  than  1.6  km  offshore  (Figure  1),  and  the  magnitude  of  the  effects  of  direct 
piling  of  water  by  the  wind  are  believed  to  be  less  than  the  range  of  error  in  tide 
measurements.  In  addition,  monthly  anomalies  of  zonal  (east  or  west)  wind  stress 
were  found  not  to  be  significantly  correlated  with  monthly  sea  level  anomalies  at 
the  5%  level  of  significance  (Table  2).  Accordingly,  the  effects  of  the  piling  up  or 
depression  of  sea  level  by  wind  stress  are  neglected  in  this  analysis. 

The  second  effect  of  wind  stress  is  that  of  sea  surface  slopes  produced  by 
offshore  or  onshore  Ekman  transport  due  to  winds  parallel  to  the  coast.  According 
to  conventional  Ekman  transport  theory,  net  transport  is  directed  90°  to  the  right 
of  the  wind  in  the  northern  hemisphere.  In  this  study,  offshore/onshore  Ekman 
transport  was  found  to  be  significantly  correlated  with  sea  level  (r  =  -0.42  in  Table 
2).  The  inverse  correlation  indicates  that  offshore  transport  results  in  decreased 
sea  level  and  onshore  transport  in  increased  sea  level.  Meridional  wind  stress  is 
also  significantly  correlated  with  sea  level  (r  =  0.43),  as  expected. 

Sea  surface  temperature  (r  =  0.61)  and  surface  salinity  (r  =  -0.35)  are  both 
significantly  correlated  with  monthly  sea-level  anomalies.  The  signs  of  the 
correlations  indicate  that  increases  in  SST  are  associated  with  increased  sea  levels 
and  increased  salinities  are  associated  with  decreased  sea  levels.  These  relation¬ 
ships  are  consistent  with  basic  considerations  of  sea  water  density  changes. 

Monthly  anomalies  of  Sverdrup  transport  were  found  not  to  be  significantly 
correlated  with  monthly  sea  level  anomalies  at  the  5%  level  of  significance. 
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B.  REGRESSION  ANALYSIS 


) 


We  have  seen  that  the  monthly  sea  level  anomalies  are  significantly  cor¬ 
related  with  atmospheric  pressure,  meridional  windstress,  offshore  Ekman  trans¬ 
port,  SST,  and  surface  salinity.  To  further  quantify  these  relationships,  a  multiple 
regression  analysis  was  performed.  The  eight  variables  input  into  the  BMDP2R 
stepwise  multiple  regression  program  (Dixon,  1975)  were  monthly  anomalies  of  sea 
level,  atmospheric  pressure,  meridional  wind  stress,  zonal  wind  stress,  offshore 
Ekman  transport,  Sverdrup  transport,  surface  salinity,  and  surface  temperature. 

The  results  of  the  regression  analysis,  presented  in  Table  3  (upper  part),  show 
that  surface  atmospheric  pressure  is  the  major  predictor  of  sea  level,  with  SST  and 
meridional  windstress  as  second  and  third  predictors.  The  remaining  variables  were 
neglible  and  their  coefficients  are  not  included  in  the  table.  Together,  these  three 

i 

variables  explain  over  68%  of  the  variance  of  the  monthly  sea  level  anomalies.  The 

i 

regression  formula  shown  in  the  table  indicates  that  the  response  of  sea  level  to 
changes  in  atmospheric  pressure  is  -1.67  em/mb  whereas  a  purely  hydrostatic 
response  would  be  -1.00  cm/mb.  This  higher  than  theoretical  pressure  response  is  in 
agreement  with  the  results  of  Saur  (1962)  and  Roden  (1960)  who  analyzed  monthly 
tide  data  at  stations  north  and  south  of  Monterey. 

Because  of  the  significant  seasonal  changes  in  the  oceanic  and  atmospheric 
regimes  near  Monterey  we  might  expect  to  observe  seasonal  changes  in  the 
processes  affecting  sea  level.  To  define  these  seasonal  changes  the  ocean  and 
atmospheric  data  contained  in  Appendix  B  were  analyzed  by  multiple  regression 
during  two  periods,  the  Davidson  Current  and  the  upwelling  seasons. 

Sea  level  changes  centered  on  the  Davidson  current  period  were  analyzed 
using  5  months  of  data  (October  through  February)  for  the  years  1963  to  1978.  As 
described  earlier,  this  is  a  period  of  weak  northerly  winds,  northward  coastal 


Table  3.  Results  of  Multiple  Regression  Analysis 


A. 

Dynamic  Height  Not 

Included 

Step 

Variable 

s  Explained 

Increase  In 

No. 

Variance 

Explained  Variance 

1 

PRES 

0.48 

0.48 

2 

SST 

0.67 

0.19 

3 

MWS 

0.68 

0.01 

SEA  LEVEL  =  0.15  -  1.67  PRES  +  2.09  SST  + 
(an)  (mb)  (°C) 


3.34  MWS 

(dynes/cnr) 


B.  Dynamic  Height  Included 


Step 

No. 

Variable 

Explained 

Variance 

Increase  In 

Explained  Variance 

1 

PRES 

0.42 

0.42 

2 

SST 

0.66 

0.24 

3 

DYNtfT 

0.72 

0.06 

4 

MWS 

0.74 

0.02 

SEA  LEVEL 
(cm) 

=  -0.32  - 

1.43  PRES  -  1.72 
(mb) 

SST  +  0.16  DYNHT  +  4.88  MWS  - 

(°C)  (dyn  cm)  (dynes/an  ) 

-  i 


current  flow,  and  frequent  cyclonic  storm  activity.  The  results  of  multiple 
regression  analysis,  shown  in  Table  4,  indicate  that  atmospheric  pressure  and  SST 
are  major  predictors  of  sea  level  during  this  period,  explaining  over  71%  of  the 
variance  of  monthly  sea  level  anomalies. 

The  second  period  analyzed  was  centered  on  the  upwelling  season  and  covered 
the  months  April  through  August  during  the  years  1964  through  1978.  The  upwelling 
season  is  a  period  of  northerly  winds,  offshore  transport  of  coastal  surface  waters, 
and  southward  California  Current  flow.  The  results  of  multiple  regression  analysis 
indicate  that  during  this  period  monthly  anomalies  of  atmospheric  pressure,  SST, 
and  meridional  windstress  account  for  58%  of  the  variability  of  monthly  sea  level 
(Table  4). 

Thus,  some  seasonal  change  in  the  processes  affecting  sea  level  is  indicated, 
with  monthly  atmospheric  pressure  and  SST  anomalies  accounting  for  most  of  the 
monthly  sea  level  variability  in  both  the  Davidson  Current  and  upwelling  seasons, 
and  meridional  windstress  explaining  an  additional  portion  of  the  sea  level 
variability  during  the  upwelling  season. 
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Table  4.  Results  of  Multiple  Regression 
Analysis  By  Season 


1 


A.  Davidson  Current  Season 


Step  Variable 

No. 


Explained  Increase  In 

Variance  Explained  Variance 


1  PRES  0.49 

2  SST  0.71 


0.49 

0.22 


SEA  LEVEL  =  0.07 
(cm) 


1.69  PRES  +  2.49  SST 
(mb)  (°C) 


Step  Variance 

No. 

1  PRES 

2  SST 

3  MWS 


Upwelling  Season 

Explained 

Variance 

0.29 

0.53 

0.58 


Increase  In 
Explained  Variance 

0.29 

0.24 

0.05 


SEA  LEVEL  =  0.58  -  1.43  PRES  +  2.27  SST  +  3.83  MWS  ? 
(cm)  (mb)  (°C)  (dynes/aiT) 


C.  SPECTRUM  ANALYSIS 


It  has  just  been  shown  that  much  of  the  variance  of  monthly  sea  level 
anomalies  can  be  explained  by  monthly  anomalies  of  atmospheric  pressure,  SST, 
and  meridional  windstress.  However,  important  variations  in  these  processes  occur 
on  time  scales  shorter  than  a  month.  To  determine  how  the  variance  of  sea  level  is 
distributed  with  frequency  over  time-periods  of  days  to  weeks,  auto  and  cross 
spectra  were  calculated  for  six-hourly  observations  of  sea  level,  atmospheric 
pressure,  and  the  meridional  component  of  wind  stress. 

Surface  atmospheric  pressure  and  meridional  wind  stress  were  calculated  as 
described  previously  on  a  six-hourly  basis  for  the  period  January  1,  1967  through 
August  31,  1976  for  a  point  approximately  14  km  west  of  the  Monterey  tide  station 
(Figure  1).  Hourly  sea  level  data  for  the  same  time  period  were  low-pass  filtered  to 
remove  the  diurnal,  semi-diurnal,  and  other  short-term  tidal  components  and  were 
sub-sampled  at  six-hourly  intervals;  a  complete  description  of  the  low-pass  filter  is 
given  by  Godin  (1966).  These  data  series  were  then  detrended  by  subtracting  their 
30-day  running  mean  to  produce  a  band-passed  series.  The  response  function  for 
the  30-day  running  mean  is  shown  in  Figure  14. 

Surface  atmospheric  pressure,  wind  stress,  and  adjusted  and  unadjusted  sea 
level  data  were  analyzed  during  the  winter  storm  period  (November  1  to  March  8) 
and  the  upwelling  period  (April  I  to  August  8)  for  the  years  1967  to  1976.  A  fast 
fourier  transform  spectrum  analysis  with  a  triangular  data  window  was  used  and 
the  seasonal  spectra  were  averaged  over  all  available  years.  The  frequency 
bandwidth  is  0.04  cpd  and  the  number  of  degrees  of  freedom  is  90  for  the  winter 
season  and  100  for  the  upwelling  season. 

The  spectral  relationships  between  sea  level  and  atmospheric  pressure  are 
shown  in  Figures  15A  and  15B  and  will  be  discussed  first.  In  the  low  frequency 
region,  the  winter  spectra  are  more  energetic  than  the  upwelling  period  spectra, 
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Figure  14.  Response  Function  for  30-day  Running  Mean  Filter 
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indicating  the  effects  of  intense  winter  storm  events.  The  largest  sea  level  and 
pressure  fluctuations  occurred  in  the  12  to  24  day  frequency  band  (on  the  event  or 
storm  time-scale).  It  is  difficult  to  say  much  about  the  spectral  peak  observed  at 
16  days.  The  analysis  scheme  used  in  this  section  was  designed  to  focus  on 
variations  with  periods  of  2  to  10  days  but  did  not  reveal  any  significant  spectral 
peaks  in  this  region.  A  recoloring  of  the  spectrum  suggests  the  possibility  of  some 
leakage  from  the  negative  filter  side-lobe  seen  in  Figure  14  into  the  frequency  band 
centered  on  0.0625  cpd. 

The  coherence  (squared)  between  sea  level  and  atmospheric  pressure  was 
found  to  be  significant  and  independent  of  frequency  in  the  upwelling  period 
(Figure  15B)  but  decreased  in  magnitude  for  periods  shorter  than  two  days  in  the 
winter  series  (Figure  15A).  The  constant  180°  phase  between  these  two  data  sets 
indicates  the  inverse  response  between  atmospheric  pressure  and  sea  level  as 
expected  from  the  hydrostatic  equation. 

To  remove  pressure  effects,  in  order  to  better  examine  the  relationship  of 
wind  stress  and  sea  level,  the  low-passed  six-hourly  sea  level  series  was  adjusted 
for  atmospheric  pressure  effects  and  detrended  in  the  manner  described  previously. 
The  six-hourly  adjusted  sea  level  and  meridional  wind  stress  series  were  then 
analyzed  and  auto  and  cross  spectra  calculated  (Figure  16A,  16B).  The  meridional 
wind  stress  also  had  a  concentration  of  energy  at  low  frequencies  with  large 
variations  occurring  in  the  12  to  24  day  frequency  band,  and  the  winter  power 
spectra  containing  more  energy  than  those  of  the  upwelling  season.  Coherence 
between  adjusted  sea  level  and  meridional  wind  stress  is  generally  low.  The  phase 
functions  provide  little  information  because  of  the  low  coherence. 
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D.  DYNAMIC  HEIGHT 


Sea  surface  temperature  and  salinity,  examined  earlier  using  multi-variable 
correlation  and  regression  techniques,  are  only  surface  samples  but  may  be 
indicative  of  the  subsurface  density  distribution.  Dynamic  height  calculations, 
however,  provide  a  direct  measure  of  the  subsurface  density  field  and  its  changes, 
and  therefore  are  a  measure  of  changes  in  ocean  circulation.  Mean  monthly 
dynamic  height  anomalies  and  long-term  means  were  calculated,  as  described 
earlier,  for  a  site  in  mid-Monterey  Bay  for  the  period  1968  through  1977.  A 
multiple  regression  analysis  was  run  which  included  monthly  dynamic  height  data 
with  that  of  the  eight  ocean  and  atmospheric  variables  described  previously.  The 
results  of  this  analysis,  shown  in  the  lower  part  of  Table  3,  indicate  that  inclusion 
of  monthly  dynamic  height  anomalies  increased  the  explained  variance  of  monthly 
sea  level  anomalies  from  69%  to  74%. 

The  relationship  between  sea  level  and  dynamic  height  was  further  examined 
in  a  seasonal  sense.  Figure  17  shows  the  mean  seasonal  cycle  of  dynamic  height  and 
adjusted  sea  level.  There  is  good  agreement  in  both  phase  and  amplitude  of  these 
curves.  The  observed  seasonal  cycle  for  dynamic  height  is  noisy  as  a  result  of 
limited  sampling  (only  5  hydrocasts  in  January  but  up  to  18  in  other  months; 
Bretschneider  and  McLain,  1979).  The  figure  shows  that  both  sea  level  and  dynamic 
height  near  Monterey  are  highest  in  winter  and  lowest  in  spring. 

Reid  and  Mantyla  (1976),  using  the  La  Jolla  tide  data  as  an  example,  showed 
that  south  of  40°  N  in  the  eastern  North  Pacific  Ocean  sea  levels  are  typically 
highest  in  late  summer  and  early  fall  and  lowest  in  late  winter  as  a  result  of  the 
annual  solar  heating  cycle.  North  of  40°  N,  however,  sea  levels  are  highest  in 
winter  and  lowest  in  summer;  this  pattern  cannot  be  explained  by  the  steric 
response  to  seasonal  heating  and  cooling.  Using  Sturges'  (1974)  data  from  Neah 
Bay,  Reid  and  Mantyla  further  demonstrated  that  maximum  sea  levels  occur  in 
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Figure  17.  Seasonal  Cycle  of  Sea  Level  and  Dynamic  Height 

(Sea  Level  Data  for  Period  1963-1978  and  Dynamic  Height 
Data  for  Period  1968-1977.  Range  of  Nfonthly  Sea  Levels 
are  Vertical  Bars) 


winter  when  inshore  northward  flow  is  strongest  and  minimum  sea  levels  occur 
during  the  southward  flow  of  summer,  thus  relating  seasonal  changes  in  sea  level  to 
geostrophically  balanced  flow.  Monterey  lies  between  these  two  stations  and  has  a 
seasonal  cycle  that  is  intermediate  between  these  regimes. 

Sea  level  and  dynamic  height  are  also  in  good  agreement  in  a  time-series 
sense.  Figure  18  shows  the  Monterey  time  series  of  weekly  mean  sea  level, 
calculated  from  the  hourly  data,  and  individual  dynamic  height  calculations 
relative  to  100,  300,  and  500  m.  The  figure  shows  that  both  sea  levels  and  dynamic 
heights  were  higher  than  normal  during  1969-1970,  1972-1973,  and  1976,  which  were 
periods  of  strong  El  Nino  activity  in  the  Eastern  Tropical  Pacific.  Sea  levels  and 
dynamic  heights  were  normal  or  near  normal  during  non-El  Nino  periods.  Because 
of  the  close  agreement  between  sea  level  and  dynamic  height,  and  the  high 
correlation  of  sea  level  at  Monterey  with  that  at  adjacent  stations,  dynamic  height 
and  sea  level  variations  both  must  be  related  to  variations  in  the  geostrophic 
current  flow  (Bretschneider  and  McLain,  1979). 
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Figure  18.  Time-Series  of  Weekly  Mean  Sea  Level  and  Dynamic  Height 
(Mean  Annual  Cycles  are  Shown  as  Dotted  Lines) 


V.  SUMMARY 


Variability  of  sea  level  at  Monterey,  California  was  analyzed  by  various 
methods  on  a  monthly,  weekly,  and  hourly  basis,  and  the  ocean  and  atmospheric 
processes  causing  this  variability  are  discussed. 

Analysis  of  13  years  of  hourly  sea  levels  indicates  that  non-tidal  sea  level 
variations  are  small  compared  to  the  normal  tide  range  in  the  area.  The  largest 
deviation  of  observed  from  predicted  hourly  sea  levels  was  39.6  cm.  A  seasonal 
change  revealed  by  monthly  frequency  distributions  of  hourly  non-tidal  sea  level 
variations  was  found,  with  observed  sea  levels  being  generally  less  than  the 
predicted  from  March  thru  May  and  greater  than  the  predicted  from  July  thru 
January. 

Monthly  sea  level  anomalies  at  Monterey  are  correlated  with  anomalies 
recorded  at  tide  stations  from  Prince  Rupert,  Canada  to  Callao,  Peru  but  are  most 
closely  related  to  events  affecting  sea  levels  in  the  group  of  stations  from 
Crescent  City,  California  to  Quepos,  Costa  Rica.  Processes  producing  the  El  Nino 
phenomenon  along  the  coast  of  Peru  also  apparently  affect  sea  level  at  Monterey. 

Multiple  regression  analysis  indicates  that  monthly  anomalies  of  atmospheric 
pressure  and  sea  surface  temperature  account  for  most  of  the  Monterey  monthly 
sea  level  variability  during  both  the  Davidson  Current  and  upwelling  seasons.  The 
meridional  component  of  wind  stress  accounts  for  an  additional  portion  of  sea  level 
variability  during  the  upwelling  season. 

Analysis  of  six-hourly  sea  level  and  atmospheric  pressure  observations  show 
that  the  winter  spectra  are  more  energetic  than  those  of  the  upwelling  season,  and 
that  most  of  the  energy  occurs  in  periods  of  12  to  24  days.  Coherence  between  sea 
level  and  atmospheric  pressure  is  significant  and  independent  of  frequency.  This 
and  the  constant  180°  phase  relationship  between  these  six-hourly  data  sets 
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reflects  the  inverse  response  between  sea  level  and  atmospheric  pressure  expected 
from  the  hydrostatic  relationship.  The  power  spectra  for  six-hourly  meridional 
wind  stress  also  show  a  concentration  of  energy  in  low  frequencies  and  are  most 
energetic  in  winter;  however,  coherence  between  the  local  wind  stress  and  sea 
level  is  generally  low. 

There  is  good  agreement  between  the  behavior  of  sea  level  and  dynamic 
height  in  both  a  time-series  and  seasonal  sense.  The  close  agreement  between  sea 
level  and  dynamic  height,  and  the  high  correlation  of  sea  level  at  Monterey  with 
that  at  adjacent  tide  stations,  suggests  that  sea  level  and  dynamic  height  both  must 
be  related  to  variations  in  the  geostrophic  current  flow. 
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APPENDIX  A 
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MISSING  HOURLY  SEA  LEVEL  DATA 

The  dates  and  times  of  missing  Monterey  hourly  sea  level  observations  are 


I 


listed  below.  The  data  series  began  July  21,  1963  and  ended  August  31,  1976. 
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APPENDIX  B 


MONTHLY  MEAN  OCEANIC  AND  ATMOSPHERIC  OBSERVATIONS 

This  appendix  presents  graphical  plots  of  monthly  means  and  monthly  mean 
anomalies  of  various  oceanic  and  atmospheric  observations  for  the  period  1960  to 
1978.  Anomalies  were  calculated  as  the  difference  between  a  monthly  mean  and 
the  long  term  mean  (1963-1978)  for  the  same  month.  Monthly  means  are  shown  as 
heavy  lines  and  monthly  anomalies  as  light  lines.  The  data  are  presented  in  the 
following  sequence: 

1)  Sea  level  (cm) 

2)  Adjusted  sea  level  (cm) 

3)  Surface  atmospheric  pressure  (mb) 

2 

4)  Meridional  wind  stress  (dynes/cm  ;  positive  northward) 

2 

5)  Zonal  wind  stress  (dynes/cm  ;  positive  eastward) 

6)  Offshore/Onshore  Ekman  transport  (metric  tons/sec  per  100  m  of 
coastline;  positive  offshore) 

7)  Sverdrup  Transport  (metric  tons/  sec  per  km;  positive  northward) 

8)  Surface  salinity  (0/00) 

9)  Sea  surface  temperature  (°  C) 
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